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There exists dark energy models that predict the occurrence of "little rip" . At the point of little 
rip the Hubble rate and its cosmic time derivative approach infinity, which is quite similar to the 
big rip singularity except that the former happens at infinite future while the latter at a finite 
cosmic time; both events happen in the future and at high energies. In the case of the big rip, a 
combination of ultra-violet and infra-red effects can smooth its doomsday. We therefore wonder 
if the little rip can also be smoothed in a similar way. We address the ultra-violet and infra-red 
effects in general relativity through a brane-world model with a Gauss-Bonnet term in the bulk and 
an induced gravity term on the brane. We find that the little rip is transformed in this case into 
a sudden singularity, or a "big brake". Even though the big brake is smoother than the little rip 
in that the Hubble rate is finite at the event, the trade-off is that it takes place sooner, at a finite 
cosmic time. In our estimate, the big brake would happen at roughly 1300Gyr. 

PACS numbers: 98.80.Jk, 04.20.Jb, 04.20.Dw 
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I. INTRODUCTION 

Recent observational evidences show that the expan- 
sion of our universe is currently accelerating [H-Q . This 
discovery has opened several possibilities to explain the 
nature of the late-time acceleration. By constructing 
reasonable models that account for this phenomenon, 
the nature of the cosmic acceleration could be either at- 
tributed to an exotic dark energy component that makes 
up roughly 70 percent of the universe [J] , or described by 
a modified gravity theory such that the predictions are 
consistent with the current observational data [f| . If the 
late-time accelerated expansion is driven by dark energy, 
the simplest setup is by adding a cosmological constant 
A to Einstein equation with cold dark matter on top of 
the baryonic sector (ACDM model), which fits very well 
the observations even though the cosmological constant 
has to be extremely fine-tuned as compared with the ex- 
pected value from quantum field theory @, @] ■ 

However, there are several other candidates that can 
play the role of dark energy, which can be character- 
ized by an equation of state w. The equation of state 
w of a given matter content is defined by the ratio be- 
tween its pressure and its energy density, i.e., the equa- 
tion of state w = P I p, where P and p are the pressure 
and the energy density of the matter, respectively. From 
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current observations, the equation of state of dark en- 
ergy is very close to —1 Nevertheless, unlike the 
cosmological constant whose equation of state is strictly 
equal to —1, the equation of state of dark energy could 
deviate from —la bit and vary with time. Further- 
more, dark energy with equation of state w < — 1 has 
not been excluded observationally [H-Q- which is also 
known as phantom energy Q . If this hypothetical phan- 
tom energy exists in the late-time universe, our universe 
will enter into a super-accelerating epoch and space- 
time singularities might occur in the future. Within a 
Friedmann-Lemaitre-Robertson- Walker (FLRW) metric, 
dark energy related future space-time singularities can 
be classified as follows depending on different types of 
divergence [§-[l2j]: 

• Big rip singularity: the singularity happens at a 
finite cosmic time t where the scale factor a, the 
Hubble parameter H and its cosmic time derivative 
H diverge [j. [15H22I ]. i.e., when t -> t s , a -> oo, 

H — > oo and H — > oo; 

• Sudden singularity: the singularity happens at a 
finite cosmic time with a finite scale factor, where 
the Hubble parameter remains finite but its cosmic 
time derivative diverges [1, [H,!!!!, i.e., when t — » 
t s , a — > a s , H — > H s and H — > oo: 

• Big freeze singularity: the singularity happens at a 
finite cosmic time with a finite scale factor where 
the Hubble parameter and its cosmic time deriva- 
tive diverge |9l. [2^ . [25H27| . i.e., when t — > t s , a — > a s , 

H — > oo and H — > oo; 
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• Type IV singularity: the singularity happens at a 
finite cosmic time with finite scale factor where the 
Hubble parameter and its cosmic time derivative 
remain finite, but higher cosmic time derivatives of 
the Hubble parameter still diverge [9j-fi~l|; 

where t s , a s and H s arc finite constants. However, in ad- 
dition to the possible space-time singularities mentioned 
above and connected with dark energy, there exist other 
possible abrupt events that are smoother than those enu- 
merated above and are also connected with dark energy: 

• Little rip event: there is no future space-time sin- 
gularity; however, the scale factor , the Hubble pa- 
rameter and its cosmic time derivative all diverge 
at an infinite cosmic time (9I l25l. I281433T ] . i.e., when 
t — > 00, a — > 00, H — > 00 and H — > 00; 

• u>-singularity: the equation of state w diverges at a 
finite cosmic time with a finite scale factor, where 
the energy density p and the pressure P vanish (3~il . 



1351 ]. i.e., when t 
w — > 00: 



t„ 



a s , p — > 0, P — > and 



• Pseudo-rip event: the Hubble parameter ap- 
proaches a constant as the cosmic time goes to in- 
finity with an infinite scale factor, where the dis- 
integration of bound structures may or may not 
occur in this case (36[, i.e., when t — > 00, a — > 00 
and H — >• H s . 

Among these events, the little rip is quite similar to 
a big rip singularity except that the former happens at 
infinite future while the latter at a finite cosmic time. 
Such event, despite avoiding a future singularity at a fi- 
nite cosmic time, will still lead to the destruction of all 
structures in the universe. Additionally, the little rip has 
previously been analyzed under four-dimensional (4D) 
standard cosmology p, [H, [H] in a dilatonic brane- world 
model [2911 . Recently, this event has been named "little 

rip" dolaa 

If our universe is filled with phantom dark energy, a 
doomsday can happen at high energies in the far future. 
These possible future behaviors have motivated the inter- 
est in searching for resolutions where the abrupt events 
can at least be appeased, and several works on possible 
future space-time singularity avoidance have been car- 
ried out (see for example the Refs. [2^, |38| - |46| ). In the 
case of a big rip, a combination of ultra-violet (UV) and 
infra-red (IR) modifications to general relativity (GR) 
can smooth its doomsday [St}- We therefore wonder if 
the little rip event can also be smoothed in a similar 
way, since it takes place in a high-energy regime where 
UV effect is required, and in a late-time universe where 
IR effect is required. Here we follow the similar procedure 
to address UV and IR effects in GR through a spatially 
flat 4D brane embedded in a 5D anti-de-Sitter (AdSs) 
bulk that is modified by both Gauss-Bonnet (GB) and 
induced gravity (IG) terms with the matter fields con- 
fined on the brane, where the GB correction manifests 



in the high-energy regime while the IG effect becomes 
significant in the late-time universe. 

Furthermore, we focus here on the self-accelerating 
branch of this model, which contains the Dvali- 
Gabadadze-Porrati (DGP) self-accelerating solution 
modified by the GB effect (DGP-GB) [H, g|, due to 
the fact that this branch gives both UV and IR modifi- 
cations to GR. Finally, the investigation on the little rip 
event can also shed some lights on the nature of such GB 
and IG modified brane-world system. 

The outline of this paper is as follows. In Sec. II, we 
consider a spatially flat brane-world within an AdSs bulk 
with GB and IG modifications, and we review the self- 
accelerating branch obtained from this model. In Sec. Ill, 
we examine if for the same matter content that would in- 
duce a little rip event in standard GR can be smoothed in 
some manner within the framework reviewed in Sec. II. 
In Sec. IV, we analyze the late-time behavior of the DGP- 
GB model filled with phantom dark energy, and further- 
more we estimate when the sudden singularity will take 
place in this case. Finally, in Sec. V, we summarize our 
results and present our conclusions. 



II. SELF-ACCELERATING BRANCH 

The generalized Fricdmann equation of a. spatially fiat 
brane embedded in an AdSs bulk with both GB and IG 
curvature effects is given by [5(3, EH : 



Vh 2 + p? 



±r r . 



^(p + \)- lH i 



(2.1) 



where the parameter p 2 is defined by 
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(2.2) 



o:(> 0) is the GB parameter which has the dimension of 
length square, 7 is a dimensionlcss parameter controlling 
the strength of the IG term, A$(< 0) is the bulk cosmo- 
logical constant, A is the brane tension and the crossover 
scale r c is given as r c = k 2 /2k 2 . Therefore the parameter 
p? is bounded as < p 2 < I/4a (see also Ref. [52| ). 

The branch with "+" sign in Eq. (|2.ip contains the 
DGP normal branch solution (with As = 0,a = 0,A = 
and 7 = 1); while the "— " sign in Eq. (|2.1[) contains the 
DGP self-accelerating solution. Here we focus on the 
self- accelerating branch with "— " sign that gives the UV 
and IR modifications to GR. To solve the Eq. (|2.I[) ana- 
lytically, it is more convenient to introduce the following 
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pi and p2 


Solutions for X 


< b < 1/4 


p2 < < pi 


Xi = -1/3 [2V1 -3bcosh(r?/3) - lj ; < pi < p 
Xi = -1/3 [2V1 - 36cos (9/3) - l] ; < p < pi 
X 2 = 1/3 [2-v/l -36cos((tt + 6»)/3) + l] ; < p < pi 
X 3 = 1/3 [2Vl -36cos((tt- 6)/3) + l] ; < p < pi 
Limiting case: 

Xi = -1/3 [2V1 - 36 - 1]; < pi = p 

X 2 = X 3 = 1/3 [Vl - 3b + 1] ; < pi = p 

Xi = X 2 = -1/3 [Vl - 3b - 1] | ft=0 = 0; = p 2 = p 

X 3 = 1/3 [2Vl - 3b + 1] | fe=0 = 1; = p 2 = p 




/Ill ,-1-1 II 
/-'J \ pi u 


Xi = -1/3 [2V1 - 3bcosh (77/3)] ; < p 
Limiting case: 

Xi = -1/3 [2V1 - 3b - l]| fc=1/4 = 0; = pi = p 

X 2 = X 3 = 1/3 [VI - 3b + 1] \ b=1/4 = 1/2; = pi = p 


b = 1/3 


pi = p 2 = -1/27 


Xi = -1/3 


(l + 27p) 1/3 - l| 


1/3 < b 


pi, P2: complex conjugates 


Xi = -1/3 


2V3b- lsinh (£/3)J 



TABLE I. Exact solutions of the modified Fricdmann equation (|2.6[) . These solutions are divided into four different regions 
depending on the value of the parameter b, and Xi, X2, and X3 correspond to three possible branches in different ranges of 
the parameter b. Here we restrict to the physical solutions where the total dimensionless energy density p is positive. 



dimcnsionless parameters |37J, |47J, |5J-|5f 

3 -fr c 

6 (l_4a^ 2 ), 
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27 7 J r2 



(2.3) 
(2.4) 

(2.5) 



Then the Fricdmann equation (|2.ip can be rewritten on 
a simpler form: 



X 3 - X 2 + bX + p = 0. 



(2.6) 



Notice that we restrict to "— " sign of Eq. (|2.ip . and X is 
positive [cf. Eq. (|2.3[) ]. Therefore only the positive solu- 
tions of Eq. (|2.6|) arc physically meaningful. 

The number of real solutions of the cubic equation (12. 6[) 
depends on the sign of the discriminant N [571 ] : 



N = Q 3 + R 2 , 



(2.7) 



where 




(2. 



It is more convenient to factorize the discriminant N as 

1 



N 
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R 2 = ^(p- pi)(p- P2), 



(2.9) 



where 



Pi 

P2 



3 1 9 



-3 \ b -9 



1 + V(l - 3&) 3 



1 - V(l - 3&) 3 



(2.10) 
(2.11) 



If AT is positive, there is only one real solution. If TV is 
negative, all the solutions are real, and if N vanishes, all 
the solutions are real and at least two of them are equal. 

After doing a careful analysis of the cubic Friedmann 
equation (|2.6[) . the analytical solutions are summarized 
in Table U (see also Fig. [1]), where the parameters 77, 9 
and £ are defined as follows: 



cosh 77 = 



R 
R 



cos( 



sinh£ = 



^Q 3 

R 



sinh i] 



sm( 



cosh£ = 



IQ 3 + R 2 



Q 3 



IQ 3 + R 2 



Q 3 



(2.12) 



(2.13) 



which give the constraints < r\ as well as < £. In 
addition, since the dimensionless energy density p is al- 
ways positive [see Eq. (|2.5[) ]. the parameter 9 is bounded 
as < 9 < 6\na X , where 



arccos 



-2 + 96 



2V(l-36) 3 



(2.15) 



In general there are three possible evolutionary 
branches for the modified Fricdmann equation (|2.6I) as 
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FIG. 1. Branches of the modified Friedmann equation (|2.6|l . 
We plot X, essentially the Hubble rate [see Eg. (|2.3[) ]. versus 
the dimensionless energy density p. For this plot, we set b = 
1/8 as an example. Notice that only positive X is allowed 
[see Eq, (|2.3p ], and the solution X2 corresponds to the self- 
accelerating branch. 



shown on Tableland Fig.[TJ but one of them is a negative 
normal branch denoted by X\, which is not physically 
meaningful since X is always positive [see the definition 
given in Eq. (|2.3[) ]. Among these positive solutions, the 
branch Xi corresponds to the self-accelerating branch: 



v 2 = i 

3 



2-s/l — 36 cos 





+1 







(2.16) 



which contains the DGP self-accelerating branch when 
the parameters A5, a, and A vanish and 7=1, which 
gives both UV and IR modifications to GR. Moreover, 
this self-accelerating branch X2 could serve to smooth the 
little rip event in standard GR to a sudden singularity, in 
the sense that the Hubble parameter remains finite when 
the universe reaches the singularity as we will show in 
the next section. 



III. SUDDEN SINGULARITY 

The late-time acceleration of the universe can be driven 
by a dark energy component with an equation of state 
w < — L i.e., a phantom energy is possible observation- 
ally [1H3|. if our universe is filled with phantom dark 
energy, this substance would lead to a super-acceleration 
phase, cau sing pe rhaps space-time singularities in the fu- 
ture |a. ll5H27l[58t . Some possible space-time singularities 
as well as smoother abrupt events within a FLRW uni- 
verse filled with phantom dark energy are enumerated in 

scc.m 

One of these abrupt events is the little rip that is quite 
similar to a big rip singularity as the scale factor, the 
Hubble parameter, and its cosmic time derivative all ap- 
proach infinity as in the big rip albeit at infinite future. 
Even though for a little rip there is no space-time sin- 
gularity occurrence at a finite cosmic time, all bound 
structures in the universe will be destroyed on its way to 
the little rip [3(J. Within a relativistic model, i.e., a 4D 



GR setup, a perfect fluid whose equation of state reads: 

P=-p-V, (3-1) 

where A is a positive constant, can induce a future little 
rip event Ip . r25l l28l . I30T ] . This kind of event was first 
found in [9| , l25l . l28j and also on a dilatonic brane- world 
model [23]. For the latter case, the author described the 
event as a big rip singularity postponed to an infinite 
future cosmic time. Recently, this kind of event has been 
named "little rip" [U. 

Instead of a relativistic GR model, here we consider 
the self-accelerating branch X2 of the modified Fried- 
mann equation (|2.6j) . because it involves UV and IR ef- 
fects which could be applied to smooth the little rip or 
at least to appease it. We remind at this regard that 
the little rip happens at very high energies where we ex- 
pect some UV corrections to GR, and at the same time 
in the far future where IR modifications to GR are also 
to be expected. This is the reason why we consider a 
dark energy component; more precisely, a perfect fluid 
described by Eq. (|3.1[) , filling the self-accelerating brane 
on top of dark matter. We will show that a future sudden 
singularity appears when the dimensionless energy den- 
sity p approaches the finite constant pi . Consequently, a 
sudden singularity will replace the little rip. 

Given that we are interested in the late-time evolution 
of the brane and the possible avoidance of the little rip, 
we assume that the brane total energy density is made 
up of non-relativistic matter as well as phantom energy 
with equation of state (|3.1I) : 



Pd = 



PmO 
„3 



PdO 



3.4 



In a + 1 



(3-2) 



where p m is the energy density of the non-relativistic 
matter, pd is the conserved energy density of the phan- 
tom dark energy (|3.1j) . p m Q and pdo are the current en- 
ergy density for matter and phantom dark energy, respec- 
tively, and the current scale factor ao is fixed as do = 1- 
In order to analyze the asymptotic future behavior of 
the self-accelerating branch X2, we need to analyze the 
future behavior of the Hubble rate given by Eq. (|2.16[) 
and its cosmic time derivative. The cosmic time deriva- 
tive of the Hubble parameter can be obtained either from 
the time derivative of the branch X2 (|2.16j) or from the 
time derivative of the cubic Friedmann equation 
Either ways lead us to the result: 



H = 



3k| 



X 2 ■ (p + P) 



7 (3X 2 - 1)2 - (1 - 36) 



(3.3) 



where we have already used the energy conservation con- 
dition for the total matter content confined on the brane: 



p + 3H(p + P) = 0. 



(3.4) 



Therefore, taking into account the total energy density 
given by Eq. (|3.2j) . the cosmic time derivative of the Hub- 
ble parameter in the self-accelerating branch is given by 



H = 



3^4 



X 2 ■ (pm ~ Ap d ) 



7 (3A 2 - l) 2 - (1 - 36) 



(3.5) 
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Notice that the denominator of Eq. (|3.5[) is always neg- 
ative in the self-accelerating branch solution, which can 
be checked using the explicit self-accelerating solution 
X2 (|2.16[) where the range of the parameter 9 is given in 
Eq. ([2.15p . Therefore, the cosmic time derivative of the 
Hubble parameter will become positive once the phantom 
dark energy dominates over the non-relativistic matter, 
more precisely when p m < Apd- 

When the dimcnsionless energy density p tends to p\, 
the self-accelerating solution X2 approaches the constant 
value (see Table Band Fig.HJ: 



Vl-36+1 



(3.6) 



consequently, the Hubble parameter reaches the constant 
value: 



H 2 



1 



36+1 -/x 



(3.7) 



From Eas. p.5|) - (|3.7p we find that the cosmic time deriva- 
tive of the Hubble parameter blows up when p — > p\, 
while the Hubble parameter, the energy density and pres- 
sure on the brane remain finite at this point. Besides, 
the cosmic time ellapsed since the present time to, till 
the singularity ocurrcnce at t a is given by the integral 
t s — to = J 1 S da/aH, which is finite since aH is a non- 
vanishing analytical function of the scale factor a between 
the two limits a s and oq = 1, i.e., a sudden singularity 
takes place at a s . Consequently, a phantom dark en- 
ergy perfect fluid (|3.1[) that induces a little rip event in 
GR is replaced by a sudden singularity, or a "big brake" , 
through the self-accelerating branch of the brane-world 
with GB and IG effects. Notice that this event is different 
from the one found in Refs. I59i 1601. 



IV. PHENOMENOLOGY AND CONSTRAINTS 
OF THE DGP-GB MODEL WITH "LITTLE RIP" 
PHANTOM MATTER 



The theoretical framework 



event. Then, the modified Friedmann equation (|4.ip can 
be further expressed as 



E(a) 2 =fi m oa 3 + 



In, 



2 £A V"do 
l + n a E(a) 2 ] E(a), 



(4.2) 



where E(a) = H/Hq, and the dimcnsionless density pa- 
rameters are defined as follows: 



_ K 4 PdO 



3H 2 



Ar^H 2 ' 



il a — —aH n 



(4.3) 



It can be easily checked that the modified Friedmann 
equation (|4.2j) evaluated at present time ao = 1 gives a 
constraint for these density parameters: 



2 v /a~(i + ft„) = 1. 



(4.4) 



Furthermore, the dimcnsionless parameters defined in the 
Eqs. (|2.3[ )- (|2.5[) can be rewritten in terms of these density 
parameters (with A5 — 0, A = and 7=1) as: 



8 a , 

X = H = 2JaTn a E(a), 

3 r c 



b=--a =4f2 a Q ro , 



8 a 
32 a 2 k% 



(4.5) 
(4.6) 



=4n r n 2 a 



n 



3A 



lna + 1 



(4.7) 



The cosmological evolution of the brane would be 
equivalent to that of a relativistic model described by 
GR, i.e., 



E(a) 2 = n m0 a 3 + n eff , 



(4.8) 



with an effective dark energy component fl c s defined as 
follows: 



In this section, we consider a spatially flat DGP brane 
embedded in a Minkowski bulk with a GB modification 
(DGP-GB), i.e., described by the modified Friedmann 
equation (|2.1[) with a vanishing bulk cosmological con- 
stant and brane tension, and 7=1. The Hubble rate of 
the branch containing the DGP self-accelerating branch 
in this case can be rewritten as 



H 



1 



-P 



-aH 2 H. 



(4.1) 



The current expansion of the universe is dominated by 
non-relativistic matter and dark energy which we model 
through Eq. (|3.2[) . the later we refer to it as "little rip" 
phantom matter as in GR it would lead to this kind of 



n>4 



3H 2 



peS = ^dO 



V3 K4A 



2 HQ\fTldo 
2y / 0~ [1 + n a E(a) 2 ] E{a). 



: In , 



(4.9) 



The effective dark energy part consists of the phantom 
dark energy as well as the higher-dimensional brane- 
world effect as shown in Eq. ([4.9[) . Furthermore, combin- 
ing the conservation equation of the effective dark energy: 



Pen + 3i?Pofr(l + w cfi ) = 0, 



(4.10) 



with the effective energy density definition (|4.9[) . we ob- 
tain the effective equation of state w e g: 



Wcff 



1 d In peg 
3 dx 



(4.11) 
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(a) (b) 



FIG. 2. The effective equation of state (|4. 11 [) of the effective dark energy (|4.9[) against the scalar factor, more precisely x = In a. 
In figure (a), we have fixed the parameter K4A = 2.75 x 10 _3 Gyr _1 , and changed the density parameter Q rc as shown on the 
plot. In figure (b), we have fixed the density parameter Q rc = e -5,56 , and changed the parameter A as shown on the plot. We 
see that the effective dark energy for the model given in Eq. (|4.2|l will induce a mimicry of a crossing of the phantom divide, 
i.e., w cross the value — 1. 




(a) (b) 



FIG. 3. The deceleration parameter q versus the scale factor, more precisely x = In a. In figure (a), we have fixed the parameter 
K4A — 2.75 x 10 -3 Gyr _1 , and changed the density parameter Q, r< . as shown on the plot. As an example, in figure (b), we have 
fixed the density parameter Q ra = e~ 5 ' 56 , and changed the parameter A as shown on the plot. 



where x = In a. This effective equation of state describes 
the evolution property of the effective dark energy p c s, 
in particular at late-time. 

In addition, the deceleration parameter q is also an 
important cosmologies! parameter that could be used to 
characterize the late-time behavior of the universe and 
to constrain the model. The deceleration parameter is 
defined by q = —'da/a , which can be rewritten as 



1 dE(a) 
E(a) dx 



(4.12) 



where the self-accelerating branch E(a) is given by 
Eqs.g3]) and (|2~T6)l . 

Before calculating the effective equation of state w s 
and the deceleration parameter q numerically, we high- 
light that both t(i c ff and q approach negative infinite val- 
ues when the universe is close to the sudden singularity 
p — > pi, since the cosmic time derivative of the Hubble 
parameter blows up at the singularity; while the Hubble 
parameter and the total energy density remain finite. 



B. Observational constraints of the model 

We assume that this self-accelerating DGP-GB model 
(|4.2[) behaves much like the standard ACDM model at 
present and on the near future. Thus we expect the den- 
sity parameters fl rc and f2 a due to the brane-world ef- 
fect to be tiny compared to the other density parameters. 
Therefore, we could pick out the results derived from the 
ACDM model to obtain some preliminary constrains of 
our model; We assume also that the density parameter 
for matter nowadays £l m o to be independent of the cos- 
mological model, i.e., we fix the parameter Cl m o = 0.27 
for the non-relativistic matter on the brane. Further- 
more, we constrain the model by using the deceleration 
parameter qo of the ACDM model at the present time: 
<7o = — 1 + 3/2 f2 m o = 0.595, which is valid for GR. 

We can calculate the deceleration parameter using 
Eqs. (|4.12[) and (|4.2[) . which at present reads 



Qo = 



2if (l 



(4.13) 
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FIG. 4. In figure (a), it shows the value of the scale factor at the sudden singularity a- s i ng 



lna s 



versus the parameter A for 



different values of the density parameter S7 rc . Here we have fixed the density parameter for matter Q m o = 0.27; while in figure 
(b), it shows the cosmic time interval before reaching the sudden singularity versus the parameter A for different values of the 
density parameter Q Tc - 



The previous relation gives a constraint among the model 
parameters, which reduces the number of free parameters 
to two (e.g. Q rc and A). In addition, since we expect 
the model to mimic ACDM cosmology reasonably well 
at the present time, the best-fit range for the parameter 
A is considered as obtained in the Ref. [30j : —2.74 x 
10~ 3 Gyr _1 < k 4 A < 9.67 x lO^Gyr" 1 . Despite that 
negative values of A arc allowed obscrvationally, we stick 
to positive values of A, such that the little rip event takes 
place in GR for a FLRW universe filled with the "little 
rip" phantom matter given in Eq. (|3.ip . 



After imposing the constraints on these parameters, we 
then evaluate the effective equation of state w e e (|4.1ip 
as a function of the scale factor numerically. Notice that 
the parameter space is picked out such that the density 
parameter f2 rc and Q a due to the brane-world effect are 
very small, and therefore the model behaves almost like 
the ACDM model at the present time. Fig. [5] shows that 
the behavior of the effective equation of state w e s with re- 
spect to the scale factor. We can see that the effects from 
the DGP-GB brane-world with phantom energy (|3.1|) will 
induce a mimicry of a crossing of the phantom divide at 
the late-time universe (see Fig. Furthermore, it is 
more sensitive to the parameter A than to the parameter 



In addition, we also calculate and plot the decelera- 
tion parameter q [see Eq. (|4.12[l ] numerically as shown in 
Fig. [3J This figure shows that our universe has entered in 
an accelerating expansion phase (q < 0) recently, roughly 
since the scale factor reached the value a ~ 0.6. Besides, 
the four curves shown in Fig. I3bl corresponding to differ- 
ent values of A, shows that the deceleration parameter 
is only sensitive to different values of A in the near past 
when a ~ 0.14. 



C. Onset of sudden singularity 

As we have seen in the previous section the little 
rip event in GR can be replaced by a sudden singular- 
ity through the self-accelerating branch of the DGP-GB 
brane-world. We can interpret this result as a smoothing 
of the little rip in the sense that the Hubble parameter 
H is finite, even though the sudden singularity happens 
before the little rip at a finite cosmic time. We further 
estimate when the sudden singularity will happen in the 
future within this DGP-GB brane-world, with a cosmo- 
logical expansion dictated by Eq. (|4.2|) and a matter con- 
tent including a phantom dark energy given in Eq. (|3.1|) . 

The sudden singularity takes place when the dinren- 
sionless energy density p — > p\ [see Fig. Q] and Eq. (|3.5l) ]. 
Therefore, in order to obtain when the sudden singular- 
ity is reached, we can directly equate the dimensionless 
energy parameter p (|4.7p to the value of the singular- 
ity point p\ (|2.10[) . and the resulting equation gives the 
scale factor x s i ng = In a s i ng where the sudden singularity 
happens. Such a scale factor, x s i ng = lna s i ng , is shown in 
Fig. [4a] for various parameter A with different density pa- 
rameters fl rc . We see as well in Fig.0a|that the larger the 
parameter A the sooner the sudden singularity happens. 
This outcome is reasonable since a larger parameter A 
generates a smaller effective equation of state w e g, i.e., a 
more negative one, as shown in Fig. I2bl which will cause 
a sooner doomsday. 

Likewise, we can calculate as well the cosmic time re- 
maining from now till the doomsday. This time scale 
can be obtained by computing the following integration 
numerically: 

f Xsins dx 

(t sing - t )H = — (4.14) 

Jo h \ a ) 

where x = In a and the self-accelerating branch E(a) is 
given by the Eqs. (|4.5[ ) and (|2.16|) . The numerical results 
are shown in Fig. I4bl where we can compare our results 
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with the age of the universe since the age of our universe 
at present is roughly of the order of the Hubble time Hq 1 . 
Thus as we can see in Fig. I4b| the sudden singularity 
would take place roughly in the far future at a cosmic 
time of the order of a hundred times the Hubble time, 
i.e., roughly 1300Gyr. 

V. CONCLUSION 

The little rip event in a relativistic FLRW universe 
described by GR can avoid the future space-time singu- 
larity at a finite cosmic time; however, it will ultimately 
lead to the dissolution of all bound structures in the uni- 
verse. Therefore, it is of interest to investigate if this 
abrupt event could be smoothed or replaced in some way. 
Since this kind of doomsday, which is induced by phan- 
tom dark energy, as the one given by Eq. p.l[) . will take 
place at high energy, where UV effect could be required, 
and in the future, where IR effect could be as well re- 
quired, we suppose that the combination of UV as well 
as the IR modifications to GR could give a resolution to 
this abrupt event. In this respect, we consider here our 
universe to be described by the self-accelerating branch 
(|2.16[) of a generalized branc-world system that is mod- 
ified by GB modification in the bulk and IG effects on 
the brane, where the GB/IG effects become significant at 
the high/late-time regime, respectively. This framework 
may provide us some insights on the avoidance of abrupt 
events in the future evolution of the universe. 

With the same kind of phantom fluid (|3.ip confined 
on the brane, we discover that a sudden singularity will 
take place in the future instead of a little rip event. Thus 
this finding shows that the little rip event is smoothed 
through a brane-world system with GB and IG effects 
as we expected, in the sense that the Hubble parameter 
remains finite during the whole evolution of the brane; 
nevertheless, the trade-off is that the sudden singularity 
will happen sooner than the little rip (but still in the far 
future), where the cosmic time derivative of the Hubble 
rate will diverge. Based on this analysis, we then study 



the late-time behavior of the brane and estimate when 
the sudden singularity will take place in the future within 
a self-accelerating branch of a DGP-GB model, under the 
assumption that the self-accelerating DGP-GB brane is 
filled with dark matter and phantom energy as the one 
given in Eq. (|3.1[) . We assume as well that the model do 
not deviate too much from the ACDM model nowadays. 
Consequently, we conclude that the sudden singularity 
will occur in a future time roughly of the order of a hun- 
dred times the age of the universe, i.e., ~1300Gyr. More- 
over, a smaller parameter A in Eq. (|3.ip will result in a 
larger effective equation of state w e g for the effective dark 
energy p c g (|4.9|) . which will push the sudden singularity 
to a farther future (see Figs. [2] and SJ. 

Finally, we point out, as demonstrated by our analy- 
sis, that in Ref. one cannot remove nor appease the 
occurrence of the little rip singularity within a DGP-GB 
setup; this is not because of the singularity itself but be- 
cause of the matter that induces it. In that case, the 
singularity was caused by a holographic Ricci dark en- 
ergy component [371 ] . 
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